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CALCULATING  GAS  FLOH  IN  A  HYPERSONIC  NOZZLE  HITR  CONSIDERATION  OF  THE 
EFFECT  OF  VISCOSITY  (DIRECT  PROBLEM) 


A.  P.  Byrkin  and  I.  I.  Mezhirov 


summary 

This  report  discusses  a  procedure  for  the  approximate  and 
precise  nuaerical  solution  of  the  problem  of  gas  flow  in  a  given 
hypersonic  nozzle  with  consideration  of  the  effect  of  viscosity  (in 
the  vicinity  of  the  boundary  layer).  The  results  of  the  computer 
calculation  of  the  flow  cf  helium  in  a  conical  nozzle  are  given. 


The  following  main  problems  make  up  the  design  complex  of  a 


hypersonic  nozzle: 
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a)  the  construction  (with  consideration  of  viscosity)  of  the 
nozzle  contour  which  provides  the  assigned  flew  in  an  inviscid 
isentropic  core  (direct  problem).  In  the  vicinity  of  the  boundary 
layer,  the  problem  is  reduced  to  increasing  the  through  cross 
sections  of  the  nozzle  calculated  ior  inviscid  flow  by  the  value 
corresponding  to  the  depth  of  displacement  of  the  boundary  layer  6*; 

b)  the  calculation  of  gas  flow  in  a  nozzle  with  a  given 
configuration  at  a  given  value  of  the  Reynolds  number  and  temperature 
of  the  nozzle  wall  (direct  problem).  The  solution  of  this  problem  is 
more  complex,  since  the  dimensions  of  the  isentropic  core  and  the 
values  of  the  Mach  numbers  in  it  are  determined  by  the  interaction  of 
the  inviscid  flow  with  the  boundary  layer,  whose  parameters,  in  turn, 
depend  on  the  characteristics  of  the  inviscid  flew. 

Among  the  studies  concerned  with  solving  the  direct  problem  for 
a  nozzle,  we  will  point  cut  reports  [1]  and  [2],  which  describe  a 
method  of  calculating  gas  flow  in  a  thin  hypersonic  conical  nozzle. 
The  integral  relationship  of  the  pulses  is  used  to  calculate  the 
boundary  layer,  the  velocity  profile  in  the  boundary  layer  is  assumed 
to  be  linear,  and  the  calculations  are  made  for  a  heat- insulated 
nozzle  wall  and  a  Prandtl  number  of  one.  A  linear  law  of  the 


k 
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dependence  of  the  aziauthal  velocity  component  on  the  polar  angle  is 
used  in  the  inviscid  core,  which  aakes  it  possible  to  consider  the 
nonunifornity  of  the  flach  numbers  induced  by  the  boundary  layer  in 
the  transverse  cross  section  of  the  nozzle. 

The  aethod  of  successive  approxiaations  is  soaetiaes  used  to 
solve  the  direct  problea  at  aodecate  supersonic  velocities  in  the 
nozzle,  whereupon  the  condition  corresponding  to  flow  of  an  inviscid 
gas  in  the  nozzle  is  used  as  the  zero  approxiaation  for  calculating 
the  boundary  layer.  The  aetnod  of  successive  approxiaations  is  also 
used  for  calculating  external  flows,  when  the  boundary  layer 
interacts  with  an  inviscid  rlow,  including  at  hypersonic  velocities 
(e.g. ,  see  [3]).  The  process  of  successive  approxiaations  usually 
converges. 

However,  one  should  oear  in  aind  that  with  a  thick  boundary 
layer  in  the  internal  preblea,  the  errors  in  deteraining  the 
thickness  of  the  boundary  layer  at  hypersonic  velocities  cause  large 
deviations  of  the  values  of  the  gas-dynamic  parameters  froa  the 
actual  paraaeters  due  tc  the  liaitation  of  the  flow  by  the  walls, 
which  aay  cause  the  iteration  process  to  be  divergent. 

This  is  illustrated  in  /ig.  1,  which  shows  the  depth  of 
diaplaceaent  of  the  boundary  layer  of  the  zero  approxiaation  60* 
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calculated  for  the  conditions  of  flow  of  an  inviscid  gas  (heliua)  in 
a  hypersonic  nozzle  with  an  apex  half-angle  of  6°,  a  ratio  of  the 
throat  radius  to  the  exit  radius  of  0.0181  (the  Mach  nuaber  of  the 
fictitious  flow  in  the  absence  of  a  boundary  layer  on  the  nozzle 
walls  llf  =  36.5),  Reu  *  848*10*  and  84.8*10*  and  with  a 

ity  i 

heat- insulated  wall.  Here  Reo i**  — — — — ,  p»0  ;  p0  and  p0  are  the 

t*o 

density  and  dynaaic  viscosity  of  the  isentrcpically  braked  gas, 
respectively;  Wm,i  is  the  aaxiauo  gas  velocity;  and  L  is  the 
length  of  the  nozzle.  In  Fig.  1  and  below, 

7_  =  ^*,  r.  =»  — ,  v  =  — .  /-  is  the  current  radius  of  the  nozzle 

r*  '* 

cross  section,  r,  is  the  radius  of  the  isentrcpic  core, 

(r,  —  rm  —  «•),  x  is  the  distance  iron  the  throat  along  the  nozzle 

axis,  and  rm  is  the  radius  of  the  nozzle  throat. 
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SO 

2S 

0 

Fig.  1. 

It  is  evident  froii  Fig.  1  mat  in  this  case,  the  calculation 
leads  to  an  absurd  result  -  the  value  of  60*  inside  the  nozzle  turns 
out  to  be  equal  to  the  nozzle  radius.  It  is  clear  that  the  iteration 
process  will  not  converge  at  auca  smaller  values  of  60*,  either. 

Belov  ve  vill  describe  an  effective  method  of  the  approximate 
and  precise  numerical  sclution  of  the  direct  problem  for  a  hypersonic 
nozzle  based  on  certain  properties  of  a  laminar  boundary  layer 
established  by  solving  the  inverse  problem  for  a  number  of  nozzles. 


1.  Figure  2  shows  dependences  of  the  dimensionless  depth  of 
displacement  of  the  laiinar  boundary  layer  \  Rtn,  on  the  number 

for  a  given  isentrcpic  contour  of  the  nozzle,  calculated  by  A.  P. 
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Byrkin  and  Yu.  N.  Pavlovskiy  for  a  nuaber  of  profiled  axisyametrical 
nozzles  with  a  heat- insulating  wall.  The  curves  which  correspond  to 
values  of  the  flach  number  of  M  *  14.9,  12.9  and  12.2  characterize 
flow  in  nozzles  with  a  tend  in  the  generatrix,  and  the  rest  -  flow  in 
nozzles  with  a  conical  section. 
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Fig. 


The  curves  were  obtained  by  the  numerical  integration  of  the 
boundary  layer  equations  by  A.  A.  Dorodnitsyn* s  method  of  integral 
relationships  (in  the  third  and  fxtth  approximations)  [4]  and  by  the 
method  of  finite  differences  [ 5 j.  The  calculations  were  conducted  for 
air  (adiabatic  index  *  *  1.4  and  a  Prandtl  number  of  Pr  =  0.75,  the 
dependence  of  the  viscosity  coefficient  cn  temperature  was  taken  from 
the  Sutherland  formula,  T0  =  10Q0°K)  ,  and  for  helium  (x  =  1.667, 

Pr  =  0.68,  the  exponent  in  tne  viscosity  law  n  -  0.647).  The  effect 


1 
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of  the  transverse  curvature  of  the  nozzle  contour  on  the  boundary 
layer  characteristics  was  not  considered  in  the  calculations,  since 
in  practice  it  can  be  disregarded  with  an  errcr  on  the  order  of 
1-2o/o  at  ordinary  Re  numbers.  This  is  indicated  by  the  results  of 
numerical  calculations  ffej,  as  well  as  data  of  special  calculations 
made  by  the  authors. 

It  is  evident  from  Fig.  2.  that  with  a  given  working  gas,  the 
dependences  obtained  for  the  eight  different  nozzles  differ  little 
from  each  other.  With  an  error  of  around  ±.10 o/o,  we  can  consider  the 
value  of  '  I  ■<<■,,,  in  hypersonic  axisy mmetrical  nozzles  with  a 
heat-insulated  wall  to  te  a  function  of  only  the  Mach  number  on  the 
edge  of  the  boundary  layer,  auu  that  it  does  not  depend  on  the  nature 
of  the  distribution  of  the  Mach  numbers  over  the  length  of  the 
nozzle.  We  will  point  cut  that  the  existence  of  this  universal 
dependence  follows  from  the  laws  of  similarity  established  by  Yu.  L. 
Zhilin  for  gas  flow  in  thin  affine-like  hypersonic  nozzles  [7]. 

£  *■ 

Figure  3  shows  the  dependence  (Tj,  plotted  for  all  eight 

''i 

nozzles.  Here  5t*  is  the  depth  of  displacement  at  a  wall  temperature 
of  'i\„  which  corresponds  to  the  case  of  heat  insulation;  6*  is  the 
depth  of  displacement  at  a  wall  temperature  of  7'a;/a,=i  is  the 
temperature  factor.  It  is  evident  from  Fig.  3  that  the  value  of 


6*/6,*  is  essentially  a  function  of  only  the  temperature  factor  Tu, 


9 
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and  it  depends  veakly  on  M*  and  the  longitudinal  pressure  gradient 
in  the  nozzle,  as  veil  as  on  the  physical  characteristics  of  the 
working  gas,  at  sufficiently  large  Mach  numbers. 


A 
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It  follows  from  the  data  given  that  in  the  first  approximation, 
the  depth  of  displacement  on  the  nczzle  wall  can  be  determined  from 
the  formula 

where  /(AUl—  Rea.,  A  (f „,)=*.  are  universal  functions  of  the  Mach 

x  £, 

number  and  the  temperature  factor  Tu„  determined  for  air  and  helium 
by  the  curves  in  Figures  2  and  i. 


2,  The  data  in  the  preceding  section  make  it  possible  to  obtain 
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a  simple  approximate  solution  to  tne  problem  of  the  gas  flow  in  a 
nozzle  with  a  fixed  configuration  during  univariate  flow  in  an 
isentropic  core  and  a  laminar  boundary  layer  on  the  nozzle  walls. 
Plow  which  is  nearly  univariate  occurs  in  conical  hypersonic  nozzles 
with  small  opening  angles  (10-15°).  The  calculation  of  flow  in  the 
core  in  the  univariate  approximation  is  often  also  sufficient  for 
profiled  nozzles  operating  in  off-design  conditions,  since  it  gives 
us  an  idea  of  the  deviations  of  tne  Hach  numbers  from  the  calculated 
values. 

We  will  write  the  equation  of  the  flow  rate  for  gas  flow  in  a 
core,  assuming,  like  in  the  calculations  discussed  earlier,  that 
there  is  no  boundary  layer  in  the  nozzle  throat: 

t-\q(M,)=F0.  (2) 

Here  F ,  and  Ft  are  the  area  of  the  isentropic  core  and  the  nozzle 
throat,  respectively;  M,  is  the  Hach  number  in  the  core;  and  q(H)  is 
the  derived  flow  rate: 


We  have 


q  (M)~~ 


«  +  i 


\ 


u  ( M 


8*  \ 
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whence,  using  (1),  we  finally  obtain 


I 


(A  ) 


\ 


X 

Re* 


/<*»,)*  </;  1 


(4> 


where  Re 


W 

"  111 


Equation  (4)  can  be  easily  solved  graphically  for  fl*  by  using 
Figures  2  and  3  for  a  nczzle  of  a  given  shape  at  the  assigned  values 
of  the  temperature  facte r  Tu  and  Re*  .  The  area  of  the  isentrcpic 
core  is  calculated  from  the  values  of  M„  (a)  using  expression  (2)  . 

It  follows  from  formula  (4)  that  the  number  M„  will  be  constant 
if  all  of  the  abscissas  of  the  points  of  contour  x  vary  in 
proportion  to  Re*  as  parameter  Be*  changes,  while  the  ordinates  r„, 
corresponding  to  them  remain  unchanged.  This  law  of  similarity 
follows  strictly  from  the  equation  of  motion  cf  a  viscous  gas  when 
the  static  pressure  is  constant  in  the  channel  cross  section  (see 
[8]). 


Ve  will  point  out  that  in  the  common  case  of  a  conical  nozzle, 
when  rw^-di\-b  •  where  a  and  b  are  constants,  we  can  obtain  the 


i 
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dependence  .r(M„)  in  explicit  fora,  here  equation  (4)  is  reduced  to  a 
quadratic  equation  in  1  x  §  and  we  obtain: 

-  (/(M.Wj  ,  >  (M,)V **7fj  ”~i' 

2a  |  R»‘„  »  4«5Re.  a 


For  a  turbulent  or  transitional  boundary  layer,  the  experiaental 
data  on  the  depth  of  displacement  can  also  be  generalized  by  a 
foraula  of  the  form 

«*  *  r» 


(e.g.,  see  [9]).  The  exponent  j  obtained  is  equal  to  0.2-0. 3.  When 
expressions  (5)  and  (3)  are  used,  we  obtain  the  following  equation 
for  determining  the  nacn  numbers  in  the  isentropic  core  of  a  nozzle 
of  a  given  shape: 


I  </(Mv) 


r,.  (x) 


•f  <M„  r.f) 


I  >1 


It  follows  froa  formula  (6)  that,  like  above,  the  nuaber  Al„  in 

_I 

the  nozzle  does  not  vary  if  the  parameters  x  and  rH,  reaain 

constant;  one  dependence  between  these  parameters  defines  a  whole 
family  of  nozzles  at  different  values  of  Re. . 


He  can  get  an  idea  of  the  precision  of  this  method  by  comparing 
the  results  obtained  using  it  with  the  data  of  the  precise  numerical 
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calculation  for  a  profiled  nozzle  (Figures  4  and  5  -  the  solid  curves 
show  the  results  of  the  precise  calculation,  and  the  broken  ones  - 
approximate).  The  contour  of  the  nozzle  vail  is  forsed  by  solving  the 
inverse  probles  by  adding  the  depth  of  displacement  to  the  contour 
calculated  vithout  considering  viscosity  (this  perfect  contour 
rj.\)  is  the  same  in  Figures  4  and  5).  The  nozzle  is  designed  to 
obtain  a  uniform  flow  uith  fl  *  ti  m  the  characteristic  exit  rhombus. 
The  wall  contour  ru  ui  la  Fig.  4  corresponds  to  the  case  of  a 
heat-insulated  vail  ( Tu  *1 ),  and  in  Fig.  5  -  to  the  value  of  the 
temperature  factor  1\  *  Q.z5.  He.  *  0. 49*10*.  The  figures  shov  the 
distributions  of  the  Hack  number  on  an  inviscid  contour  obtained  by 
the  precise  numerical  calculation.  They  also  shov  the  curves  of 
r,  (,vi  and  M.  (vi,  obtained  ay  solving  equation  (4)  at  given 
ru-  Hi-,,  Tu  •  It  is  evident  that  tne  approximate  dependences  differ 
little  from  the  precise  cnes  over  the  entire  length  of  the  nozzle, 
although  the  depth  of  displacement  of  the  boundary  layer  at  the 
nozzle  exit  even  exceeds  the  radius  of  the  inviscid  core. 


i 


DOC  =  0868 


PAGE  16 


approximation  results  in  a  clow  pattern  close  to  the  true  one.  As  the 
preceding  section  shows,  we  can  use  the  data  obtained  by  solving  the 
approximate  equation  (4)  as  the  zero  approximation. 

The  results  of  numerical  computer  calculations  of  the  flow  of 
heliun  in  a  conical  nozzle,  the  aain  characteristics  of  which  are 
given  at  the  beginning  cf  the  article,  at  Ri*0/  848*10*,  530*10*, 

84.8*10*  and  two  temperature  conditions  on  the  wall  -  TW^TU„  and 
7^  =  7,,  (T0  is  the  stagnation  temperature  of  the  gas)  are  given 

below. 


The  boundary  layer  was  calculated  by  the  method  of  generalized 
integral  relationships  in  the  third  approximation. 

Axisymmetrical  gas  flow  in  an  inviscid  isentropic  core  whose 
contour  was  determined  by  the  equation 

rj  'c)  rv(x\—  (.r).  (7) 

was  calculated  by  the  method  of  cnaracteristics  using  a  specially 
written  program  for  solving  the  direct  problem.  There  were  at  least 
125  points  on  each  characteristic.  It  was  assumed  that  there  is  no 
boundary  layer  in  the  nczzle  throat.  In  order  to  avoid  the  detailed 
calculation  of  the  transonic  section  of  the  nozzle,  it  was  assumed 
that  radial  gas  flow  at  fl  *  1.01  occurs  immediately  after  the  throat. 
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The  Mach  nuaber  on  the  boundary  of  the  isent topic  core,  obtained 
froa  foraulae  (5)  and  (2)  in  ua  zero  approxiaation,  was  used  to 
calculate  the  depth  of  displaceaent  of  the  boundary  layer  of  the 
first  approxiaation.  Then  new  values  of  the  radius  of  the  inviscid 
core,  etc.,  were  deterained  until  the  approxiaation  yielded  virtually 
the  saae  results  with  identical  satisfaction  cf  relationship  (7). 

Since  the  iterations  liuctuate  around  an  unknown  liaiting  value 
at  supersonic  velocities  in  an  inviscid  core  (this  follows  froa  the 
aain  gas-dynaaic  relationships  -  the  area  of  the  channel  increases  as 
the  S  nuaber  increases  -  ana  the  tact  that  the  value  of  ^  )•  Re«r  is 
an  increasing  function  ct  the  Hach  nuaber) ,  "daafing“  was  used  to 
iaprove  convergence:  the  radius  or  the  inviscid  core  in  the  i-th 
approxiaation,  which  is  used  to  calculate  the  flow  of  the  inviscid 
gas,  was  calculated  froa  the  foraula 

''ip.o,  •  'V,), 

where  the  daaping  coefficient  {  was  considered  to  be  equal  to 
0.5-0.25. 


The  process  always  turned  out  to  be  convergent,  and  the  nuaber 


of  approxiaatlons  required  did  not  exceed  four.  The  data 
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corresponding  to  a  heat-insulated  wall  were  used  as  the  zero 
approxiaation  for  the  nuaerical  solution  of  the  problea  with  the 
boundary  condition  Tw—T„  .  Two  approximations  were  necessary  in 
this  case. 

Figure  6  shows  the  distribution  of  Hach  nuabers  over  the  radius 
of  the  nozzle  exit  secticn  obtained  by  these  calculations,  it  is 
evident  that  at  Reo/.  *848«106  and  530«106,  the  Hach  nuabers  in  the 
inviscid  core  markedly  decrease  with  distance  from  the  nozzle  axis. 
Hhen  Re<>i  *  84.8«10*,  the  flow  in  the  core  is  clcse  to 
unidiaensional. 


i 

» 


L 
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The  calculated  and  experimental  distributions  of  Hach  numbers  on 
the  nozzle  axis  at  Re-./  *  U4.BOO*  are  compared  in  Fig.  7  (the 

experimental  study  was  conducted  by  V.  Ya.  Eezmenov  and  1.  I. 
ftezhirov)  *. 


Footnote:  ‘The  Knudsen  cumber  calculated  for  the  nozzle  radius  did 
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not  exceed  (1-2)  «10~2  is  the  experiments,  which  indicates  the 
validity  of  the  comparison  of  the  experimental  and  calculated  data. 
End  footnote. 


This  figure  also  shows  the  curve  M, <Vi  ,  which  corresponds  tc  flow 
of  an  inviscid  gas  in  a  conical  nozzle.  The  agreement  of  the 
calculated  and  experimental  data  is  satisfactcry.  The  effect  of 
viscosity  on  the  flow  in  the  nozzle  is  characterized  by  the 
difference  in  the  actual  values  of  the  Nach  numbers  from  the  values 
of  ,M corresponding  tc  radial  flow  of  an  inviscid  gas. 


Figure  8  shows  the  velocity  profile  and  the  stagnation 
temperature  profile  in  the  boundary  layer  in  the  nozzle  exit  section 
for  Rent  =  530«106  and  a  heat- insulated  wall  (the  variable  , 


p  y 

,, - n  r*  V  Re,  - 

^max  P06  f , 


is  proportional  to  the  distance  from  the  wall  y;  and  Tu,.  are  the 
values  of  the  velocity  and  stagnation  temperature  on  the  outer  edge 
of  the  boundary  layer,  respectively).  The  graph  shows  the  value  of 
corresponding  to  the  depth  of  displacement.  It  is  evident  that 


the  depth  of  displacement  differs  insignificantly  from  the  thickness 
of  the  boundary  layer. 
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In  conclusion,  we  will  point  cut  that  with  consideration  of 
viscosity,  the  procedure  given  here  for  the  numerical  calculation  of 
gas  flow  in  a  nozzle  with  a  given  shape  is  valid  for  an  arbitrary 
nozzle  with  a  sufficiently  smooth  contour  and  an  arbitrary  gas  at 
arbitrary  boundary  conditions  on  the  wall.  It  is  only  limited  by  the 
requirement  of  the  validity  or  using  the  boundary  layer  equation 
(i.e. ,  for  example,  the  condition  of  the  absence  of  shock  waves 
interacting  with  the  boundary  layer  in  the  nozzle,  the  condition  of 
sufficient  smallness  of  the  effects  of  the  rarefaction  of  the  gas). 
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